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Abstract 

We calculate the cross sections for the production of Higgs particles in association 
with a photon in e + e~ collisions, e + e~ — ► 7+Higgs, allowing for the longitudinal 
polarization of the initial electron and positron beams. We consider the associated 
production of both the Standard Model Higgs boson, and the neutral CP-even 
and CP-odd Higgs particles of its minimal supersymmetric extension. Complete 
and compact analytical expressions are given, and the size of the cross sections is 
illustrated for energies which will be reached at future e + e~ colliders. 
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1. Introduction 



The search for Higgs particles and the study of the electroweak symmetry breaking mech- 
anism |IJ is one of the most important goals of future high-energy colliders. Once Higgs 
particles are found, it will be of utmost importance to make a detailed study of their 
basic properties. In particular, the measurement of the couplings of the Higgs bosons to 
the other fundamental particles will be a crucial test of the nature of the Higgs bosons. 
In this respect, future e + e~ colliders || will play a major role: the clean environment in 
these machines will allow a rather precise determination of these couplings. 

In e + e~ collisions, the couplings of the Standard Model (SM) Higgs boson H° to the 
massive gauge bosons W and Z can be precisely measured in the main production chan- 
nels: the bremsstrahlung process e + e~ — > H°Z || and the WW/ZZ fusion mechanisms 
e + e~ — > W*W*/Z*Z* — > H°v e v e /H°e + e~ {§]. The Higgs couplings to fermions can also 
be measured in some mass range: for Higgs masses below the WW threshold, the strength 
of the Higgs couplings to r leptons and charm quarks, relative to the dominant coupling 
to bottom quarks can be determined ||. The large Higgs Yukawa coupling to heavy top 
quarks can be measured in the production process e + e~ — > itH° || or if the Higgs is heavy 
enough to decay into top quarks, in the process e + e" — > H°Z — > ttZ 0. Finally, the 
trilinear Higgs self-coupling can be accessed via the double Higgs production processes 
e+e- -> ZH°H° and e+e" -> W*W* -> v e v e H°H° |]. 

Another set of important Higgs boson couplings consists of the H°gg, H 0/ -f-f and H°Z / ~f 
vertices. These couplings do not occur at the tree level but are induced by loops of heavy 
particles. Because of the Higgs interaction being proportional to the particle masses, the 
particles do not decouple for very large masses, and these vertices could therefore serve to 
count the number of particles which couple to the Higgs boson. The H°gg vertex |J can 
be accessed through the decay H° — > gg in e + e~ collisions for Higgs bosons in the mass 
range below ~ 140 GeV, or in the fusion reaction gg — > H° which is the main production 
mechanism of Higgs particles at LHC. The if °77 and H^Z^i couplings can be accessed in 
the decays H° — > 77 \W, |TT|] and H° — > Z7 [T2|] which typically have branching ratios 



of O(l0~ 3 ) and could be measured at the LHC. The if°77 width can also be determined 
directly by means of the 77 fusion process [with the photons generated by Compton back- 
scattering of laser light] at high-energy e + e~ colliders ||13|| , while the H°Z / ~f coupling can 
be measured in the rare decay of the neutral vector boson Z — ► H 0/ ~f JO], [14], [L5| , if H° is 
light enough. 

A precise determination of these couplings could help to distinguish between the SM 
Higgs boson and Higgs particles predicted by some of its extensions such a two-Higgs 
Doublet Model (2HDM), supersymmetric theories (SUSY) [p]] or some other type of New 
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Physics ||16|| . For instance, in the Minimal Supersymmetric extension of the Standard 
Model (MSSM), the Higgs sector is enlarged to contain two doublets of scalar fields 
leading to a spectrum of five Higgs particles: two CP-even Higgs bosons h and H, a CP- 
odd Higgs boson A and two charged Higgs particles H ± . In the decoupling regime [ O , 
where the Higgs bosons H, A and H ± are very heavy, the lightest MSSM Higgs particle h 
has exactly the properties of the SM Higgs boson, except that its mass is restricted to be 
smaller than M h < 140 GeV |l9f| . If the SUSY particles are too heavy to be produced 
directly at the colliders, the only way to discriminate between the SM Higgs boson H° 
and h in the decoupling limit is to look at the Higgs couplings to Z7 and 77. While in the 
SM these couplings are induced by heavy fermion and W boson loops only, the MSSM 
provides additional contributions from loops involving charginos and sfermions [^0[ plfl . 



An alternative way to have access to these induced Higgs^y7 and Higgs- Z7 couplings 
is to consider the process 



e e 



#°7 • (1-1) 



In the SM, this process proceeds through s-channel 7*7./?° and Z*^H° vertex diagrams, 
but additional t-channel vertex and box diagrams involving iy/neutrino and Z/electron 
exchange also occur (Fig. la). This process has been discussed some time ago in Ref. p2| 
and more recently in Ref. for unpolarized electron and positron beams. Since it is a 
higher-order process in the electroweak couplingf], the cross section is rather small, but 



the signal is very clean allowing for a resonable hope to isolate these events [^5] if enough 
luminosity is collected at a future high-energy collider. 

In this paper, we rederive the cross section for the process eq. (1.1) in the Standard 
Model, but allowing for the longitudinal polarization of the initial e + and e~ beams. 
In addition, we consider the associated production of a photon and a Higgs boson in 
supersymmetric theories. In the MSSM, the associated production of the CP-even neutral 
Higgs particles 

e + e~ —>■ h'j 

e+e" -> #7 (1.2) 

and the production of the CP-odd Higgs boson 

e + e" -> A7 (1.3) 

will receive additional contributions coming from SUSY loops. For the CP-even Higgs 
bosons one has additional s-channel 7, Z exchange vertex diagrams involving loops with 

1 A similar process is the loop induced double production of Higgs particles, e + e~ — > $$, which has 
been discussed in the SM and the MSSM recently [El. 
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charged Higgs bosons, charginos, squarks and sleptons as well as t-channel vertex and box 
diagrams involving chargino/sneutrino and neutralino/selectron loops (Fig. lb). For the 
pseudoscalar state A, because of CP invariance, the process is mediated only by s-channel 
vertex diagrams involving chargino and fermion loops, as well as t-channel vertex and box 
diagrams involving chargino/sneutrino and neutralino/selectron loops (Fig. lc). 

The rest of the paper is organized as follows. In the next section, we analyze the 
associated Higgs-photon production in the Standard Model. In section 3, the case of the 
neutral CP-even and CP-odd Higgs bosons of the minimal supersymmetric extension is 
discussed. Some additional technical material is given in a short Appendix. 



2. Associated production in the SM 

2.1 Decomposition into form factors 

In the Standard Model, the process e + e~ — > ^H is described to lowest order by the 
Feynman diagrams depicted in Fig. la. There are s-channel 7 and Z exchange vertex 
diagrams involving virtual W and heavy fermion loops [mainly top and bottom quark 
loops, since the Yukawa couplings to the light fermions are very small], as well as t- 
channel vertex diagrams involving neutrino and Z/electron exchanges [corrections to 
the H°ee vertex] and W/neutrino and Z/electron box diagrams. Additional contributions 
come from diagrams where the Higgs boson is emitted from the virtual Z line and with 
Z-7 mixing through fermion or W boson loops; diagrams with the mixing between the 
Z boson (or the photon) and the Higgs boson give negligible contributions. 

Allowing for the polarization of both the initial electron and positron beams, and 
summing over the polarizations of the photon, the differential cross section of the process 
can be written as 

dcos# V [ ; 2s 16tts ^ 

where y/s is the cm. energy and 9 the scattering angle of the photon. Neglecting the 
electron mass, the amplitude M. can be decomposed into the following sum of amplitudes 

M = £ 12 K c i (2-2) 

i=l,3 v=+- 

where the form factors Cf sum all diagrams of Fig. la: 

Cf = C2 ± + Cf ± + C e t ± + C } r ± . (2.3) 
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C 7 and C z denote the contributions of the 7 and Z pole vertex diagrams, C e the t-channel 
H°ee vertex corrections and C box the contributions of the box diagrams; the contributions 
due to the Z-j mixing are included in C z . Denoting by p± the momentum [all momenta 
are taking to be incoming] of the initial e ± beams, p 7 and e 7 the photon momentum and 
polarization vector, the set of standard matrix elements Af is given by 

Af = v(p+) (1 ± 75) (A PvP- - fa e ~fP~) U (P-) 
» 2 = v (p+ ) ( 1 ± 75 ) (jL, p~, .p+ - fa e 7 .p+ ) u (p_ ) 

= !7( P+ )(l± 75 )0£ y )«( P _). (2.4) 

Summing over the photon polarizations, the squared amplitude can be written as 

1 



A^ 



with 



T,\M\ 2 = 



pol 



(i6tt 2 )< 



E 



E 



£"YD rpVV 



(2.5) 



ij 



EA, + (A/) t 

pol 



*T = EA l -(A-) t = 



pol 



:i + A_)(l-A+)T i 



-(l-A_)(l + A + )T y 



T 1 — ^ = T 1 "' — = 







(2.6) 



where A_ and A + are the longitudinal polarizations of the initial e~ and e + beams, re- 
spectively. In terms of the usual Mandelstam variables s = (p_ + p+) 2 , u = (p_ + p 7 ) 2 
and i = (p+ + p-y) 2 , the matrix T reads 



T 



2s 



/ u 

2* 
V 2u 2t 4 



2«\ 

2 



(2.7) 



/ 



In fact, as we will discuss later in more details, QED gauge invariance requires that the 
form factors Cf vanish after summing over the contributions of all Feynman diagrams. 
This is obvious from eq. (2.4): while the form factors Af 2 vanish if the polarization of the 
photon e 7 is replaced by its momentum p 7 , it is not the case for Af ; the transversality of 
the photon therefore implies that it is the amplitude Cf which should vanish. Therefore, 
one can drop the contributions of A3 and the differential cross section of the process 
7 if can be written in the simple form 

A " S ~ Ml 1 f l + A_)(l-A. 



e + e 



dcos# 



64tts 



1 

+ (1-A_)(1 + A + ) 



u 



\ 2 + t 2 \c+\ 2 



„ 2 I/0-I2 1 4. 

u \U 1 I + 1 



2 'CVI 2 



(2.8) 



5 



where the scattering angle 9 is related to t by t = — (s — M| r )(l — cos#)/2. The unpolarized 
cross section is obtained by simply setting A± = 0: 



J^-unpol. 



Ml 



dcos# 64stt (16tt 2 ) 2 



u 



{\ct 



2 + \cr\ 2 )+t 2 



2^ 1 2 + 1^2 



(2.9) 



2.2 s— channel diagrams 

Let us now discuss the contributions of the various diagrams to the form factors Cf and 
start with the s-channel vertex and Z-photon mixing contributions. The contributions 
of the fermions and W bosons to the r y r jH Q or Z / ~fH° vertex diagrams can be decomposed 
into the following tensorial form [all momenta are taken as incoming and V = Z, 7]: 



(2.10) 



Neglecting the electron mass, for on-shell final photons only the form factors G\ and G3 
contribute to the amplitude. They are related to the form factors Cf 2 7 by 



c7 ± 
Cf ± 



C 9 7± 



e 1 

2 s 3 

ez ± 



AswCw s — Mf 



Gi 



(2.11) 



2s^ with s 2 ^ 



l-cfy = sin 2 6> W - 



where e is the electric charge, z + = — 1 + 2$^ and z' 

Although only the form factors Cf 2 7 contribute to the e + e~ — > if°7 cross section, we 
also give the expressions for C-f' 7 to check explicitly that QED gauge invariance is indeed 
fulfilled: 



C 7± 



G 



z ± 



e 1 
2 s 

e z 



Gl 



Mh g , 



± 



qZ _ s-Ml qZ 



(2.12) 



4s w c w s — Mf 

Note that in the previous expressions we have omitted the finite width in the Z boson 
propagator for simplicity. 

The form factors G] ,z are obtained by summing bosonic and fermionic contributions 
to the V^H° vertices: 

,2 



Gl = 



G? = 



e 3 M w 



m 



f 



e 3 M, 



w 



CW Sy/ 



,Z,W 



M w 



f Ivl w 



f pf 

1 1 



(2-13) 
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where Qf, l( and mj are the electric charge, weak isospin and mass of the fermion / [in 
practice only top and bottom quarks] respectively, and N c is the color factor, N c = 1 for 
leptons and N c = 3 for quarks. 

The explicit calculation of the form factors Gj' 3 was done in the Feynman gauge where 
the W boson contributions can be split into the g^ v and charged Goldstone parts. The 
ultraviolet divergent amplitudes have been reduced from a complicated tensorial form to 
scalar Passarino-Veltman two- and three-point functions using the package FORM 



27j]. These scalar functions are defined in the Appendix and their complete analytical 



expressions can be found in Ref. [2S| for instance; their numerical evaluation has been 



performed using the package FF |29| 



The fermion loops give the same contribution F[ 3 to the 7*7^° and Z*jH° vertices 
since the difference due to the different 7// and Zff couplings has been factorized out 
in eq. (2.13): 

F( = -fl + (2m 2 -M 2 H - s)C - 2s(2C n + C 21 ) + 2(s - M 2 H ){C 12 + C Tc 
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Fi = C + AC 12 + AC 23 (2.14) 

with the Cq and three-point functions defined as 

C y = Cii(s, 0, M 2 H ; m 2 f , m 2 f ,m 2 f ) . (2.15) 

Note that the fermionic contribution to the Z—j mixing diagrams is proportional to the 
final photon momentum squared and therefore vanishes for photons on the mass shell. 

The W boson loops involving true vertex diagrams [with the exchange of W bosons, 
charged Goldstones and ghosts] and two-point functions [involving the quartic couplings 
between two scalars and two vectors] give different contributions to the j*jH° and Z*^fH 
vertices because of the complicated structure of the trilinear gauge boson and quartic 
gauge boson scalar couplings which are different for the Z and the photon. Summing over 
all diagrams, the W boson contribution reads 

F?> w = (^ + 6\(4C 24 -B 13 )+M 2 H (-7C + C 11 + C 12 ) 



+s(5C + C u - C 12 ) - B 12 + B 



' M- 



2 



J 23 



F 3' W = 4(-f- + 6)(C 12 + C 23 ) + 16Co (2.16) 



for the 7*7if° vertex, and 
1 M 2 

F ?' W = oT^( 1 - 2c w)^C 24 -B 13 ) + Mlc 2 w (7C -C 11 -C 12 -2C 23 ) 
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+M 2 H (-C + C 11 - 
-c 2 w [s(5C + C n 



C 23 ) + 2M^(c 2 w -l)C 

C\2 + 2C*21 — 2C23) + 32C*24 - 

C21 + C23) - B12 + -B23 



£1, - 65 



13 



B 



23 



2c 2 w ) + 2(1 - 6c : 



w J 



(C 1 2 + C2 3 )+4(l-4c^)Co 



(2.17) 



for the Z*7if° vertex. The two-point and three-point functions, B^ [note that here, 
the Bij functions are different from the standard ones given in for instance] Cy- are 
defined as 



-B12 
-B13 
-B23 



B (s;M^,M^) 

B Q (M 2 H -M 2 W ,M 2 W ) 

B (0;M^,M^) 



(2.18) 



The expressions in eq. (2.17) include the W boson contribution to the Z-7 mixing dia- 
grams, which is non-zero only in the case of the form factor F^ ,w and reads 

2B 23 = 2B (0;M 2 ,,M^) . (2.19) 



" 1 I mix 



Finally, we note that while the fermionic contributions are gauge invariant by them- 
selves, i.e. that one has C 3 + = for fermions, the W contributions are not. Gauge 
invariance of the W contributions is obtained only when the t-channel vertex and the box 
diagrams are included. 



2.3 t— channel vertex and box diagrams 

The contributions of the t-channel H°e + e~ vertex corrections come from W^/neutrino 
and Z/electron loops. Only the "transverse" part of the gauge bosons contribute since 
the longitudinal or Goldstone part is proportional to the very small electron mass. The 
contribution of these diagrams to the form factors Cf can be written as 



c ' ± = 7- 



Mi 



w ^w± 



+ 



Mz 
4 c 3 



z± 



w 



+ crossed 



with the VT/neutrino contributions 



(2.20) 



AT 



C 12 (m 2 e ,t,M 2 H] M^,0,M^) 




w± 



A 



w± 







(2.21) 
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and the Z/electron contributions 

Ai± = {z^fC^ml^Ml-MlmlMl) 

Af± = Af ± = . (2.22) 

The crossed terms are obtained by substituting t — * u in the functions C\ 2 of eqs. (2.21) 
and (2.22). Because of the t-channel electron exchange, terms proportional to ln(mg) 
appear and one therefore has to keep a finite value for the electron mass in the argument 
of the C\2 functions. This dependence on ln(m e ) will disappear after adding the W and 
Z box contributions. 

Finally, the contribution of the box W and Z diagrams to the form factors Cf°* can 
be written as 

, box± _ e 4 M w r „ w+ . , n . e 4 M. 



44- 



Br + crossed] + B?± . (2.23) 



The terms from the FT/neutrino box contribution read: 

BY + = 4(D + D n -D 13 + D 23 -D 25 ) 
5f + = A(-D 12 + D 13 + D 23 -D 26 ) 

3s(-D 12 + D 13 + D 25 + D 26 - D 23 - D 24 ) 

+ 7t{-D 23 + ,026) + u{3D 13 + 7D 25 - 7D 23 ) - 20D 27 ] . (2.24) 



Br = l 



Since the W boson is left-handed, there is no contribution from the W boxes to the form 
factors ■ 



}W- 



B™- = . (2.25) 

The four-point functions D and have to be understood as 

Dij = D^ml ml M 2 H , 0, s, u; M 2 W , 0, M 2 W , M 2 W ) (2.26) 

in the expressions given above, while in the crossed contributions one has to interchange 
the photon and Higgs boson masses and interchange u <-> t in eqs. (2.24) and (2.26) 

D tj = AiK,m2,0,M|,s,t;M^,0,M^,M^) . (2.27) 

The Z/electron box contributions in eq. (2.23) are given by 

5f± = 2(z ± ) 2 (-£> n + D 12 + D 22 - D 24 ) 
B^ = 2(z ± ) 2 (D 22 - D 26 ) 

Bl± = r f [s(D 22 - D 24 + D 25 -D 26 ) + 2D 27 ] (2.28) 



where 



D tj = D l3 (m 2 e ,M 2 H ,m 2 e ,0,t, U ;M 2 ,m 2 e ,m 2 e ,M 2 ) . (2.29) 

There is no crossed contribution here since the photon does not couple to the Z boson 
directly. Again, we have kept the electron mass in the arguments of the functions to 
avoid infrared singularities; these logarithmic singularities cancel exactly, as it should be, 
those which appear in the t-channel H°ee vertex contributions. 

2.4 Results 

We have verified explicitly that the sum of the W and the Z boson contributions from the 
various diagrams is gauge invariant, i.e. that the form factor C3 in eq. (2.3) indeed vanishes 
when all contributions are added. This cancellation occurs as follows: the Z/electron 
contributions to the H°ee vertex and the box diagrams cancel each other, C3 ± + C3 ox ± = 
0. Furthermore, the W boson contributions to C] ~ + C$ ~ is zero, and as discussed 
previously the contributions to of the W boson induced H°e + e~ vertex and box 
diagrams are zero since the W boson is left-handed. Finally, the form factor C% vanishes 
after summing the W boson contributions of all diagrams. This feature, together with 
the fact that the sum of all contributions is both ultraviolet and infrared finite, provides 
a very strong check of the calculation. 

We have also compared our results with previous calculations of the e + e~ — > H°^y 
cross section in the unpolarized case |E2^ |23|. Our analytical results agree completely 



with the result s0 of Ref. |23|], provided that their fermionic contribution to the s-channel 
7*7 H° and ZjH vertices is multiplied by a factor of two. Our result for the fermionic 
contributions agrees with the one obtained in Ref. JTT], 15, 14] for the Z — > H 0/ y decay 



and also with Ref. |22| if the color factor for quarks is included in the fermionic sum. 



We have also compared numerically our results with those of Ref. [P4 , who calculated 
the e + e~ —>■ H 0/ y cross section using the Feynman gauge. Our numbers agree with those 
given in Tab. 4 of Ref. |22|] up to a few percent; the small difference is probably due to a 
slightly different choice of input parameter values. 

The cross sections for the process e + e~ — > H 0/ y are shown in Fig. 2 as a function of 
the Higgs boson mass for two center of mass energies \/s = 500 GeV and 1.5 TeV typical 
of future high-energy e + e~ colliders. At 500 GeV, the unpolarized cross section [solid 
curve] is of the order of a ~ 0.05 fb for light Higgs masses and increases slightly when 

2 The results of Ref. |^3) have been obtained using a non-linear gauge and the individual contributions 
are therefore different from ours; however the sum of all contributions is gauge invariant and can be 
compared with our final result. 
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approaching the WW threshold; it then drops out quickly with increasing M H due to the 
lack of phase space. At yfs = 1.5 TeV, the cross section for light H° drops by a factor of 
~ 5 compared to the previous case, but the decrease with increasing Mh is slower. 

Fig. 2 shows also the effect on the cross section of polarizing longitudinally the electron 
and positron beams. With left-handed longitudinally polarized electrons [dashed lines] 
the cross section a{e/ j e + — > H 0r y) is approximately two times larger, and with left-handed 
electrons and right-handed positrons [dotted lines] the cross section o{e^e\ — ► H°j) is 
approximately four times larger than in the unpolarized case. The longitudinal polariza- 
tion of both the electron and positron beams is therefore important to increase the cross 
section; for energies in the 300 GeV range the cross section can reach values of the order 
of 0.5 fb for M H ~ 100 GeV. 

Fig. 3a exhibits the dependence of the cross section on the center of mass energy 
for several values of the Higgs boson mass. The cross section increases rapidly with the 
opening of the phase space and then decreases near y/s ~ 350 GeV close to the tt threshold 
[the W and top quark contributions interfere destructively, and the top contribution is 
maximal near the ti threshold] . The cross section drops smoothly with increasing energy; 
despite of the presence of t-channel vertex and box contributions, it scales approximately 
as 1/s at high energies. For Mh ~ 100 GeV, the cross section can be maximized by 
running the collider at energies around 250 GeV. 

Finally, Fig. 3b shows the angular distribution dcr/dcos^ for a cm. energy of 500 
GeV and a Higgs boson mass Mh = 100 GeV. The distribution is forward-backward 
symmetric and does not depend very strongly on the Higgs boson mass. 

With the yearly integrated luminosity of £ ~ 100 fb -1 expected at future high-energy 
e + e _ colliders, one could collect a few tens of H 0r y events in the course of a few years, if 
longitudinal polarization is available?]. The signal, which mainly consists of a photon and 
a bb pair in the low Higgs mass range or WW and ZZ pairs for Higgs masses larger than 
~ 160 GeV, is extremely clean. The backgrounds should be very small since the photon 
must be very energetic and the bb or WW/ZZ pairs should peak at an invariant mass Mh 
[these rare events will be searched only once the Higgs boson is found in one of the main 
production channels and the Higgs mass would be precisely known by then]. Therefore, 
despite of the low rates, the clean signal gives a reasonable hope to isolate these events. 

3 Note that at the CERN collider LEP2 with a cm. energy of y/s ~ 180 GeV, the expected luminosity 
which will be available, C ~ 0.5 fb -1 , does not allow to produce any H°j event. 
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3. Associated production of MSSM Higgs bosons 



In the MSSM, the associated production of the CP-even Higgs bosons h, H and the CP- 
odd Higgs boson A with a photon receives extra contributions from loops involving charged 
Higgs bosons and genuine supersymmetric particles. More precisely, the production of the 
CP-even Higgs particles, e + e~ — > $7 with $ = h,H, proceeds through the same diagrams 
as for the SM Higgs boson substituting the Higgs couplings to W, Z bosons and fermions, 
plus those of Fig. lb: (i) s-channel 7, Z exchange vertex diagrams involving loops built 
up by charged Higgs boson, chargino, squark and slepton loops, (ii) t-channel diagrams 
correcting the <J>ee vertex involving chargino/sneutrino and neutralino/selectron loops, 
and (Hi) box diagrams involving chargino/sneutrino and neutralino/selectron states. 

For the associated production of the pseudoscalar state A, because of CP invari- 
ance, there is no contribution from W/Z bosons, charged Higgs bosons and sfermions in 
the vertex diagrams and no contribution from W/Z bosons in the box diagrams. The 
process e + e~ — > A^y is therefore mediated only by s-channel vertex diagrams involv- 
ing chargino and fermion loops, as well as t-channel vertex and box diagrams involving 
chargino/sneutrino and neutralino/selectron states (Fig. lc). 

The differential cross section for e + e~ — > $7 with $ = h,H and A is given by the 
same formulae as in the SM, eq. (2.1), with Mh replaced by M$. Again, since by virtue 
of gauge invariance the form factor Cf vanishes, the cross section with polarized e + e~ 
beams reduces as in eq. (2.8) to 

u 2 \C+\ 2 + t 2 |C+| 2 " 



+ (1-A_)(1 + A + ) 

In spite of this, we give Cf to explicitly check gauge invariance. The form factors Cf 2 
and Cf are different for the production of CP-even and CP-odd Higgs production. We 
therefore discuss these two cases separately. 

3.1. CP— even Higgs boson production 

A. s— channel vertex corrections 

The one-loop 77$ or ^7$ induced vertices can be written in the same tensorial form as 
in eq. (2.10). For the CP-even Higgs particles, the form factors Cf are given in terms 
of the amplitudes Gf as in the SM by eqs. (2.11-2.12). In the MSSM, the various loop 



d a 
d cos# 



s - M$ 



QAstt (lQn 2 ) 2 



1 + A_)(1-A + ) 



\ 2 + t 2 



(3.1) 
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contributions to the form factors Gf are given by 



Gl = 



G 



e 3 M, 



w 



e 3 M, 



Fr + F^ + -J2^Q}N C ^LF/ 
/ m w 



e 

s w 



IV 



CwS w 



Ff' W + F t Z ' H+ - £ *Qf N c - 2s wQf ) F} 

f 



w 



s w 



Ft x +J2 N cQfF; 
f 



z,f 



(3.2) 



F^ and F( denote the W boson, the charged Higgs boson, 



where F?' w , F^ H+ , 

the chargino, the sfermion and the fermion contributions, respectively. Note that while 
the contributions of the fermions, charged Higgs bosons, charginos and sfermions are 
separately gauge invariant, the contributions of the W bosons are not [as in the SM] 
and the contributions of the t-channel vertex and box diagrams are needed to insure the 
gauge invariance of this subset of diagrams. The contributions of the various loops are as 
follows: 



W boson loops : 

The contributions of the W bosons to the 7*7$ vertex are 



F[ 



g<s>G+G- — 

C W 



M W C - 4C 24 + B 13 + g 9VV M|(-6C + C n + C 12 ) 



+s(5C + C n - C 12 ) + 24C 24 - B 12 - 6B 13 + B 23 

-4-^*g+g-4 (C12 + C 23 ) + g*vv& (2C + 3(Ci2 + C23)) 
c w 



(3.3) 



and the contributions to the Z*7$ vertex are 



F{ 



z,w 



gzvv[Mlc 2 w (-6C + C u + C 12 + 2C 23 ) + M|(-C n + C 23 ) 

+2M W {1 - c 2 w )C + c 2 w (-l + s(5C + Cn - C 12 + 2C 21 - 2C 23 ) + 32C 24 

-B 12 - 6B 13 - B 23 ) - 1/2 - sC + sCn + sC 21 - sC 23 + B 12 - B 23 



+9&G+G- 



M w (l - c^)C + (1 - c^/2)(4C 24 - B 13 ) 



4 \(Ac 2 w - 1)C + (6c 2 w - 1)(C 12 + C 23 ) 



+9*G+G- 2 {c^ -2){C l2 + G 



23, 



(3.4) 



Note that the bosonic (W, G + ) two-point function contribution to photon-Z boson mix- 
ing, Fi' w \ mix = 2B 23 g<$>vv is included. The two- and three-point functions B and C are 
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the same as in eq. (2.18) with Mh — > M$. 

At first sight, the couplings of the Higgs bosons to vector bosons and charged Gold- 
stones, which are given in Tab. la, look quite different. However, in the MSSM the two 
couplings are intimately related; for instance in the case of the light CP-even Higgs bo- 
son [a similar identity holds for the heavy CP-even Higgs boson], one can writeQ these 
couplings in terms of the ratios r h = M^/M^ and r H = M^/M^ 



Qhvv = sin(/3 - a) 



r H (r H 



9hG+G~ 

and therefore 



cos 2(3 x sin(/3 + a 



9hG+G~ 

9hvv 



rn - r h )(r H + r h - 1) 



1/2 



r H + r h -l 



1/2 



x 



r h [r H 



th - r h 



1/2 



Ml 

Ml 



(3.5) 



(3.6) 



Inserting this identity, which is also valid in an arbitrary two-Higgs Doublet Model (see 
for example pPf), in eqs. (3.3-3.4), one then recovers the standard W boson contributions 
eqs. (2.16-2.17) with an overall additional factor g&vv- 

Fermion loops: 

The contributions of the fermion loops are as in the SM, modulo a global factor 
for the couplings of the Higgs bosons to fermions relative to the SM Higgs couplings (see 
Tab. lb): 

,f _ 9<s>ff - 



fi 



2 

9*ff 



1 + (2m) 
Co + 4C 12 



Ml - s)C - 2s(2C u 
-4C 23 



C 21 ) + 2(s-Mi)(C l2 + G 



2:-! i 



(3.7) 



The three-point functions C , Cy are defined as in eq. (2.15) with M H — >■ M$. 

Charged Higgs boson loops: 

The contributions of the charged Higgs bosons to the 77$ and ^7$ vertices are the 
same as those of the charged Goldstones, and read 



FT 



Fr' 



—g>S>H+H- + 2g$VV (4C 2 4 



B 



13 



1 



c 



2~9^>H+H- + 2g<$VV 4(Ci2 + C23) 



(31 



4 These expressions are valid only at the tree-level, however the ratio should not be altered by the 
radiative corrections in the Higgs sector if the latter are consistently included in both the Higgs boson 
masses, the mixing angle a and in the trilinear ghG+G- coupling. 
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F Z,H+ 



2 

W 



2g<s>vv + -^-9<!>h+h- (4C; 



-IV 



--24 



B 



13 j 



{2c 2 w - 1) 2g^ vv + —g* H +H- 4 (C 12 + C 



-■23 ) 



(3.9) 



The couplings g^n+H- are the same as the charged Goldstone couplings g^c+G- given in 
Tab. la at the tree level; however, if the radiative corrections to the Higgs boson masses 
are incorporated, one should also use the corrected coupling g<$>H+u- given (|3.51|) , to 
insure the gauge invariance of the result. The two- and three point functions are defined 
in this case as 



B 

d 
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'ij 



B (M 2 ;M 2 H± ,M 2 H± 



^{s,0,Ml;M 2 H± ,M 2 H± ,M 2 H± ) 



(3.10) 



Chargino loops: 



The chargino contributions Fj ,x with i = 1, 3 to the V^7$ vertices can be cast into 
the compact form 



F 



E /*( 

j,fc=l,2 L 

E 9 



x 



X, 



A.B=L,R 



B 



(3.11) 



The couplings of the charginos to the Higgs bosons are given in Tab. lc, while the chargino 
couplings to the photon and the Z boson read 



L/R r j L/R 



,L/R /{ 



SwCw) 



(3.12) 



The matrices 0' L I R are defined as 

0% = -VaV£-±V a V£ + 5 ij 8ly , 0[f 
where the U and V matrices which diagonalize the chargino mass matrix can be found in 



—UixUji - -U i2 U* 2 + 5ijS 2 w 



(3.13) 



Ref. i3l| . The functions fi/gi = fi/gi(mi } m 2} m 3 ) are: 



fi 
9i 

h 

9s 



mi 



m 2 
+m 3 



2m 2 m 3 C - M|(C + C n + 2C 12 + 2C 23 ) 
- s(C + 3Cn - 2C 12 + 2C 21 - 2C 23 ) - 4C 24 + 
Ml(Cn + C 12 + 2C 23 ) + s(Cu - C 12 + 2C 21 - 2C 23 ) + 8C 24 - 1 



M|(C 12 + 2C 23 ) - s(2C u - C 12 + 2C 21 - 2C 23 ) - 4C 24 + 1 



2mi C + Cu + 2C 12 + 2C 
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m 2 (Cu - C i2 ) + m 3 (Ci 2 + 2C 23 ) 
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(3.14) 



The arguments of the three-point functions C and Cj,- are specified according to 



Cij — Cij(s, 0, M$] m \+i m x+' m x + ^ ' 



(3.15) 



Sfermion loops: 

Finally, the squark/slepton contribution to the ^7$ vertices can be written as 



9&f j f k 9vfafj Si 

j,fe=l,2 



(3.16) 



This term comes only from Higgs emission from the sfermion lines since the contributions 
of the two- and one-point Z-7 mixing diagrams cancel each other. The functions Sj = 

sArrif ,m; , m ? ) read: 

' V Jj' Jfc' Jk> 



s 3 



2B (Ml; m) ,m\) - 8C 24 (s, 0, M|; m 2 h ,m\,m\) 



-8[C 12 (s,0,M|;m 



m|) + C 23 (s, 0, M|; m|, , m\)\ . (3.17) 



/fe' 7fe- 



The squark couplings to the Higgs bosons, including mixing J^J between left- and right- 
handed sfermions0 are 



9*hh 
9<s>hh 
9®hh 



ct L cos 2 e f + c^ 

C* R co S 2 9 f + Ct L sm*e f 



sm 2 9 f 1 



2C RL cos 6* j sin 9 f 



2C RL cos9j sin 9f 
C%l (cos 2 0, - sin 2 e f ) + (Ct R - Cf L ) cos 9 f sin 0, 



(3.18) 



with 



C 



LL 



M 



c w 
M z 



Z'jf 



QfSw)gf 



m 



Mi 



/ $ 
-01 



IV 



cw 
m f 



(QfS 2 w )9 



ni 



2Mi 



<.i.> ..„$" 



/ $ 
0i 



^4/04 + ^03 



(3.19) 



where Af and are the soft-SUSY breaking trilinear term and the Higgs-higgsino mass 
parameter, respectively. The factors gf are given in Tab. 2a, and in the limit of zero- 
mixing the couplings reduce to those given in Tab. 2b. The photon coupling to sfermions 

5 This mixing is proportional to the fermion mass and in practice, is non negligible only for the partners 
of the third generation fermions. 
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is just g 7 f.f. = SijQf, while the couplings of the Z boson to sfermions in the case of mixing 
are given by 



9zhh = Dll cos 2 9f + Drr sin 2 9f 
9zf 2 f 2 = D RR cos2 Of + d ll sin 2 9 f 
9zhh = ( Drr ~ Dll "> sin6 f cos e f 



with 



D 



LL 



SwCw 

B. t— channel vertex corrections 



(// - Q;S 2 W ) and Drr 



SwCw 



-QfS 



w) 



(3.20) 



(3.21) 



In addition to the contribution of the W/ v and Z/e SM like loops, there are two additional 
contributions to the t-channel $ee vertex diagrams: one with chargino/sneutrino loops 
and another one with neutralino/selektron loops. In the case where the mixing in the 
selectron sector is neglected, all these diagrams do not contribute to the form factors C{%: 



C{% = . 

There is, however, a contribution to the form factors Cf which reads 



(3.22) 



c e 3 + 



'w L 



M, 



w 



M? 



9 g*vv fl VV (M w ,0,M w ) - ^ gs, vv (z + ) 2 f{ vv (M z ,m e ,M z ) 



-w 



+ E E 9e X fv 9e X +v9t x + x - fir(-m xp m ,-m xt ) 



A=R,L j,k 

M z 
2c 



w 



E 9 ex +u 9e X +v9®vi> fl SS {mz, -m +,m ) 



3 



1 

9 

z A=R,L j,k 

Mz 



+ ^ E E ^exScx, 9e X l~e L /ff^ , m e ~ L , 771*0 ) 

+ 7— E 5e X ?ez, #e X °ez, ^ez, J 3 ( m e L , ™ x ° , "lg L ) 



+ crossed (3.23) 



cr 



Mz 
M z 



9*vv (z-) 2 fi vv (M z ,m e ,M z ) 



A=R,L j,k 



+ ^E ^exjgfl Se X °e, ^e K e« //^K^^o^eJ 



+ crossed (3.24) 
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with A = L,R when A = R,L and with the functions /j" = f$'{m\, m 2 , m 3 ) given by 



ffSS 
J3 



— (-^12 

2 



fVV 



ffs 



C±2 
1 



fx 



ffS 
SB 



m 3 C 



12 



involving the three-point functions Cy specified as 

= Cijiml, t, Ml; ml, ml, m\) 



(3.25) 



(3.26) 



for the direct diagrams; for the crossed diagrams one has to make the substitution t — > u in 
the previous equation. The couplings of the Higgs bosons to vector bosons, charginos and 
neutralinos are given in Tab. 1, and those to sleptons in Tab. 2b. The only remaining cou- 
plings to be defined are the electron-chargino-sneutrino and electron-neutralino-selectron 



couplings; normalized to gw 



1/2 



M w , they are given by 



ex- v 



n Sw AT 

2 — lVi2 , g ex Oe L 



-N i2 - — ^-iVji . 
cw 



(3.27) 



The matrices V and N which diagonalize the chargino and neutralino mass matrices can 
be found in Ref. [j3l[l. 



C. Box corrections 

The box diagrams involve the contributions form W /neutrino and Z/electron loops as in 
the case of the SM Higgs boson, plus chargino/sneutrino and neutralino/selectron loops. 
The contributions to Cj 30 ^ read: 



a 



boxH 



„3 
S W 



^ fI VVV (M w , 0, M W , M W ) - ^fl vvs (M w , 0, M W , M W )) 



V" nfSSSf \ 

f 2^ ^ 9e x°eL 9e X °e L 9$e t e L // {m SL , m x o , m~ eh , m~ eh ) 



E E 9 ex tu 9 ex tv gL+ x - f!l fS (-m xt ,m P , -m + -m +) 



A=_R,L j,k 



+ cros. 



+ - 



M 



4r 3 



Z (^ + )// /yV/ K,M z ,M z ,m e 



+ E E g 9e X p L 9£ x ° x ° 91a S (™ x ° , m e L , m SL , m x o ) 



A=R,L j,k 



M 



+ E ^7.9e X fi> 9 ex i» 9*** fi (-m xt ,mp, m p , -m xt ) 



a 



box- 



M 



IV 



"X-k 



+ cros. 



+- 



M z 
4<w 



[z-)f t ffVV (m e ,M z ,M z ,m e 



+ E E^e fl ^x°e«^ x o x o ^/i 55 (m o,me fl ,m e - K ,m x o) (3.28) 

A=R,L j,k J 

where the sums run over the two chargino and the four neutralino states. All the couplings 
have been previously defined, and the box functions f/ xyz = /j z (mf,m|,m|,m|) are 
given by 



jfSSS 

f, 



fSSS 
2 



ffSSS 
J3 

ffVVS 
Jl 

ffVVS 
J2 



fVVS 



ffVVV 
Jl 

rfVVV 
J2 



fVVV 



ffffS 
JlL 

ffffS 
JlR 

fffS 



2L 

fffS 



f 
J2R 



fffS 
3L 



ffffs 
J3R 



= -D 

= -D 
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2.3 



D 
D 



25 



26 



= -[t(D 23 -D 26 )+u(D 23 -D 25 ) + 2D 27 ] (3.29) 

= 2(D + D n ) 
= 

= i[-2w(A) + D u ) - (s + 2t)D 12 + (s + 2t + u)D 13 + (s + t + u)D 23 

+ sD 2A - (s + u)D 25 - (s + t)D 26 ) + AD 27 \ (3.30) 

= D + D u - 2D 13 + 2D 23 - 2D 25 
= 2[-D 12 + D 13 + D 23 - D 26 ] 

= \[u(D Q + D n ) - (s - t)D 12 + (s-t + u)D l3 - 2(s + 2t + 2u)D 23 

- 2sD 2A + 2(s + 2u)D 25 + 2(s + 2t)D 26 - 12D 27 ] (3.31) 

= -m^Do + D u ) + m A (D 13 - D 23 + D 25 ) 

= m 3 (-D 23 + D 25 ) 

= m A (-D 23 + D 26 ) 

= m 3 (D 12 - D 13 - D 23 + D 26 ) 

= \\i. m i ~ m 4 )sD 12 - (m 1 - m A )sD 13 + (m^s + u + t) - m A s)D 23 
+ (mi — m A )sD 2A — ((mi — m 4 )s + miu)D 25 

- ((mi - m A )s + mit)D 2 Q + (4m x - 2m A )D 27 ] 

= -m 3 (mim 4 A) - s(D 23 + D 2A - D 25 - D 26 ) - 2D 27 ) (3.32) 
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with the four-point functions 

D i:j = Aj( m e> ml, Mf , 0, s, u; m\, m\, m\, ml) . 
The remaining box functions read 

fffVV : 



fffVV 
J2 



ffvv 
ffss 



fl 

fffSS 
J2 



f: 



ffss 



2(-D u + D l2 + D 22 -D 24 ) 
2{D 22 - D 2% ) 

-m x m 4 D§ + s(D 22 - D 24 + D 25 - D 26 ) + 2D 27 
—D 22 + D 2 4 

-D 12 + D 13 - D 22 + D 26 

-[mim 4 D - s(D 22 - D 2A + D 25 - D 2% ) - 2D 27 ] 



with the four-point functions 

Dij = Dij(ml, Ml, m\, 0, t, u; m\, m\, m\, ml) , 



and 



Jfss 

9lL 

Jfss 

9lR 

n ffSS 
92L 



m 1 (D + D 11 + D 12 + D 24 ) 
m A (D X2 + D 2A ) 
m 1 (D 13 + D 26 ) 
g g ss = m 4 (D 26 ) 



Jfss 

93L 

Jfss 

93R 



m 1 (-u(D + D u + D l2 + D 24 ) - t(D 13 + D 26 ) - 2D 27 ) 



with 



= 7>m 4 (-u(D 12 + D 24 ) - tD 26 - 2D 27 ) 



Dij = Dij(m 2 e , 0, m 2 e , Ml, u, t; m\, m\, m\, ml) 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



The contribution of the crossed diagrams can be obtained by simply interchanging t «-> u 
in the relevant expressions given above. 

The gauge invariance of the complete result has been checked explicitly. The sum of 
VF/neutrino and Z/electron contributions to the t-channel vertex and box diagrams and 
the contributions of the W bosons to the s-channel vertices are gauge invariant as in the 
SM. For the neutralino/selectron and the chargino/sneutrino contributions, only the sum 
C e + C box is gauge invariant. 



20 



3.2. CP-odd Higgs production 



In the case of the CP-odd Higgs bosons, fewer diagrams contribute to the associated 
production with a photon compared to the case of the CP-even Higgs particles. Indeed, 
because of CP invariance, there is no pseudoscalar couplings to vector bosons, to charged 
Higgs bosons and to the SUSY partners of the light fermions. The mixing between 
left-handed and right-handed third generation sfermions will induce a Afif 2 coupling, 
however, CP-invariance still imposes the sum of all sfermion loops to vanish. 



A. s— channel vertex corrections 

The one-loop 77A or Z'-fA vertices can be written in the same tensorial form as in 
eq. (2.10). However, for the CP-odd Higgs boson, the form factors Cf 2 are different be- 
cause now the amplitude G$ gives a non-vanishing contribution. Using the same notation 
as previously, we have for the Z'jA vertex 

r z± ez± 1 / r z , r z\ 

° 2± = t^mi T G * ] 

<** = 4^7^^-^^), 0.39) 
and for the 77A vertex 



cy* = 


e 
2 


1 

s 


(-GJ ± Gl) 


c 2 7± = 


e 
2 


1 

s 


{-Gl =F Gl) 


= 


e 
2 


1 

s 


(g1 - s -f* 



Gl . (3.40) 



Only heavy fermion and chargino loops contribute to the form factors G] and G z 
They are separately gauge invariant and are decomposed according to 



Gf = ^ 1 E2Q/iV c ^f(/3 / -2^Q / )^ 



— (3.41) 



s w 



,3 

The fermion loop contribution to the VAj vertex F[ is given by 

F( = F( = 

Ft = -g Aff C (s,0,M 2 A ;m 2 f ,m 2 f ,m 2 f ) (3.42) 
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with the couplings g^ff in Tab. la. The chargino loops yield 



Fif + = E hA- m xp- m x + ^- m xt) + 3i^{-rn x+kJ -m + ,-m +) 



E [/«( 



)"*( 



m +, -m + , -m 4 

Afc A j Aj 



X 



fe+xr + g vxtxz) ( 9 Mx7 - 9 A 



Xk X-j 



(3.43) 



with / 13 and g 13 given by eqs. (3-14) and the new functions 

h = 2m 1 (C + C 11 ) 

g 6 = 2(m 2 Cn - m 2 Ci 2 + m 3 Ci 2 ) . 

The arguments of the C^- functions are specified as in eq. (3.15) with M$ = M A . 
B. £— channel vertex corrections 



(3.44) 



The t-channel Aee vertex corrections are built-up only by chargino/sneutrino and neu- 
tralino/selectron loops since there is no AW coupling. Again, these loops do not con- 
tribute to C\ 2 and the expression of the contribution to Cf is simpler than in the case 
of the CP-even Higgs bosons, since in the absence of slepton mixing, there is no All 
coupling: 



/-<e± 
°1,2 



cr 







(3.45) 



E E 9e X +v 9e X +v9 a Ax + x7 (~ m xt ' m *> ~ m xt ) 



a=R,L j,k 



X k 



+ \ E E ^ x o e - L 9e X »e L <fe° fH S (m x °,m iL ,m x o) 



a=R,L j,k 



s 3 



2 E E 9e X °e R 9exle R <?A X ° X ° fU (™X° > m e«, ™ x o ) 
a=R,L j,k 



+ crossed 
+ crossed 



(3.46) 



with the functions fif S [a — L when a = R and vice versa] are given by eq. (3.25) with 
the same functions as in eq. (3.26) with the replacement M$ — > Ma, and the crossed 
contributions are obtained by interchanging u and t. 
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C. Box corrections 



Only the chargino/sneutrino and neutralino/selectron boxes contribute to the production 
of the pseudoscalar. Using the same notation as in the case of the CP-even Higgs bosons, 
the C; box± amplitudes are given by 

„4 



a 



boxH 



E E 9e x +i> 9 ex U 9A x j x - Ha ' ^ '^t ' ~ m xt ) 



a=R,L j,k 



+ cr. 



e 4 1 



C. 



box- 



+ 3-^ E E 9e X p L 9e X \e L 9%£,& > > , ™ x ° ) 

e 4 1 

3" 9 E E 9ex?S« 2e X °e« ^x?x£ 9ia ( m x° > m ^ > m e R , ™ X \ ) 



(3.47) 



where sums run over the two chargino and the four neutralino states [again a = L,R 
when a = R,L}. The various couplings and the box functions /:•" and gf have been 
given previously eqs. ( |3.32| - p737D ; the crossed contribution is obtained by the interchange 
u <->• t. 

Again, we have checked explicitely that the sum of the contributions C e + C hox is 
indeed gauge invariant. 



3.3. Results 

Before discussing our numerical results, let us first shortly describe our parametrization of 
the MSSM parameter space. Besides the four masses, M h , M H , Ma and M H ±, the Higgs 
sector is described at the tree level by two additional parameters, tan j3 and a mixing angle 
a. Due to supersymmetry constraints only two of them are independent and the two inputs 
are in general taken to be tan/3 and Ma- Radiative corrections JTSL |T!|, the leading part 
of which grow as the fourth power of the top mass and logarithmically with the common 
squark mass, change significantly the relations between masses and couplings and shift 
the mass of the lightest CP-even Higgs boson upwards. These radiative corrections are 
very important and should therefore be included in the analysis. We will, however, only 
include the leading part of this correction which in the simplest case can be parametrized 
in terms of the quantity [18 

_ 3G F 

V2n 2 sin 2 /? 

The CP-even Higgs boson masses are then given in terms of the pseudoscalar mass Ma 
and tan (3 as 

1 




(3.48) 



Mi+M'+e 
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=F V( M i + M z + e ) 2 - 4MiM| cos 2 2/3 - Ae{M\ sin 2 /3 + Mj cos 2 /3) 



(3.49) 



and the mixing angle a by 



M 2 A +Ml it , 

tan 2a = tan 2/3— = £ -, <a<0. 3.50 

P Mj - M| + e/ cos2/3 ' 2 - ~ V 7 

Once tan /3 and are chosen and the leading radiative correction is included in a and 
M h /M H , all the couplings of the Higgs bosons to fermions, gauge bosons, charginos, 
neutralinos and sfermions are fixed; these couplings are given in Tables 1 and 2. There 
is, however, one exception: the trilinear Higgs couplings also receive radiative corrections 
which are not completely mapped into the corrections to the angle a. The only trilinear 
couplings which we will need in our analysis are the CP-even Higgs boson couplings to 
charged Higgs bosons and Goldstones [which are the same at the Born level]; to leading 
order they are given by []33 



nn . , n \ e cos a cos 2 /3 

9hH+H- = cos 2/3 sm(/3 + a) + —r — 

M| sin p 

e sin a cos 2 (3 



Ml sin (3 



e 



9hh+h- = -cos2/3cos(/3 + a) + 

cc 

9hg+g- — — cos2/3cos(/3 + a) + — j sin a sin/3 (3.51) 

M z 



9hG+G~ — cos 2/3 sin(/3 + a) + — j cos a sin (3 



M 



2 
Z 

e 



It is these expressions which have to be used in order to obtain the same results for the 
lightest CP-even Higgs boson h and for the SM Higgs boson H° in the decoupling limit, 
and to insure the gauge invariance of the final results. 

In the MSSM, there are two charginos xtl^ = 1, 2] and four neutralinos \i [i — 1-4]. 
Their masses and their couplings are given in terms of the Higgs-higgsino mass parameter 
fi, the ratio of the vacuum expectation values tan/3 and the wino mass parameter M2. 
The bino and gluino masses are related to M2 by Mi ~ M2/2 and rrig ~ 3.5M2, when 
the gaugino masses and the three coupling constants of SU(3)xSU(2)xU(l) are unified 
at the GUT scale. The squark and slepton masses are given in terms of the parameters 
fi and tan/3, as well as the left- and right-handed scalar masses Mi and Mj r , which 
we will take to be equal, and the soft-SUSY breaking trilinear coupling Af. In the case 
of the third generation sfermions, mixing between left- and right-handed states is taken 
into account; in pratice it is important only for top squarks. 
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The total cross sections for the associated production of the lightest CP-even Higgs 
boson, e + e~ — > /ry, in the unpolarized case are shown in Figs. 4-7. The cross sections for 
longitudinally polarized electron and positron beams will not be displayed: as in the SM 
case, they can be obtained by simply rescaling the figures by approximately factor 2 for 
left-handed electrons and a factor 4 for left-handed electrons and right-handed positrons. 

Fig. 4 shows the unpolarized cross sections at a centre of mass energy y/s = 500 GeV 
as a function of M h for two values of tan/3, tan/3 = 2 and 50. In Fig. 4a(4b), the SUSY 
parameters have been chosen in such a way that charginos, neutralinos and sleptons will 
(not) be kinematically reachable at a 500 GeV e + e~ collider. The solid lines show the 
total cross sections including the SUSY contributions, while the dashed lines include only 
the standard and Higgs contributions as it is the case for a two-Higgs doublet model 
(2HDM) where the MSSM relations have been implemented. 

For low Mh the cross sections are smaller than in the SM, especially for large tan/3 
values. This is due to the fact that the main contributions are coming from the W bosons 
loops, and in the MSSM or 2HDM, the W boson couplings are suppressed by a factor 
sin(/3 — a) compared to the SM Higgs coupling. The suppression can be rather drastic, 
especially for large values of tan/3, where the cross sections drop to very small values. 
This can be best seen in Fig. 5 where the cross sections are plotted against tan /3 for two 
values of Ma = 50 and 500 GeV. With increasing M h) the cross sections increase and in 
the 2HDM they reach, as it should be, the SM cross section values in the decoupling limit 
where M h is maximal [in practice for M A ~ 500 GeV which leads to M h ~ 90 GeV for 
tan/3 = 2 and M h ~ 120 GeV for tan/3 = 50]. 

The contributions of the SUSY particles interfere destructively with the contributions 
of the standard particles in most of the parameter space. The only exception is when the 
standard contributions give a very small cross section, as in the case of very high tan/3 
values far from the decoupling limit. In the scenario of Fig. 4a and Fig. 5a, where SUSY 
particles can be found directly at a cm. energy of 500 GeV, the difference between the 
MSSM and the 2HDM is very large, more than a factor of 2. Unfortunately, since the 
cross section is smaller than in the SM, the signal will be harder to isolate compared to 
the case of the standard Higgs boson. For the scenario where the SUSY particles are too 
heavy to be produced directly at the chosen cm. energy, the SUSY contributions are very 
small and they will be very difficult to be detected. 

This is best illustrated in Fig. 6, where the relative difference between the SM and 
the MSSM cases in the decoupling limit Ma 3> M z . In this limit h and the standard 
H° have practically the same couplings and it is very hard to distinguish between the 
SM and the MSSM. This difference is plotted against M 2 in Fig. 6a and \x in Fig. 6b, for 
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tan/3 = 2,50 and for several values of fi and M 2 , respectively. For small values of the 
parameters M 2 or /i, charginos are light enough to give sizeable contributions, especially 
near the production threshold. For larger values of these parameters, the SUSY particles 
are too heavy and since their couplings are not proportional to their masses, they decouple 
from the production amplitudes rendering the difference between the MSSM and the SM 
contributions very small. 

In Fig. 7, the cross section e + e~ — > h'y is shown as a function of the cm. energy for 
tan /3 = 2, 50 and Ma = 50 and 500 GeV. As expected from the SM case, the cross sections 
decrease smoothly [except when some threshold for the real production of a particle is 
crossed] with increasing ^/s. Therefore, to optimize the cross section, one needs to operate 
the collider at energies around ^fs = 200 GeV. This is about the LEP2 cm. energy; but 
at LEP2, the luminosity is so small that no event can be expected. 

Finally, Figs. 8 and 9 display the cross sections for the associated production of the 
heavy CP-even Higgs boson H [Fig. 8] and the pseudoscalar Higgs particle A [Fig. 9] at 
y/s = 500 GeV as a function of the Higgs masses. The same scenarios as in Fig. 4 have 
been chosen. Except for light H and A and for small values of tan (3, the cross sections 
are very small. Again, the difference between the MSSM and 2HDM is large only when 
the SUSY particles are within the reach of the collider. 

4. Conclusions 

We have calculated the cross sections for the production of Higgs particles in association 
with a photon in e + e~ collisions, allowing for the longitudinal polarization of the initial 
e~ and e + beams. We have considered the case of the Standard Model Higgs boson, 
e + e~ — > 'yH , and the case of the neutral CP-even and CP-odd Higgs particles of the 
Minimal Supersymmetric extension of the Standard Model, e + e~ — > 7/1, jH and 7A 
We have given complete and compact analytical expressions, and made careful checks of 
the gauge invariance of the results, as well as a comparison with earlier results for the 
associated production of the SM Higgs boson in the case of unpolarized beams. 

We have then illustrated the size of the cross sections for energies which will be reached 
at future e + e~ colliders. In the SM, the cross sections are in general small but the signals 
are rather clean. With an integrated luminosity of £ ~ 100 fb _1 expected at future high- 
energy e + e~ colliders and with the polarization of both electron and positron beams, 
which increase the cross sections by a factor of 4 compared to the unpolarized case, one 
could hope to isolate the signals despite of the low rates. 

The cross sections for the associated production of the lightest MSSM CP-even Higgs 
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boson are in general smaller than for the SM Higgs boson except in the decoupling limit 
where they reach the SM values. The contribution of the SUSY particles are significant 
only when these particles are light enough to be produced directly at the collider. For 
SUSY particles beyond the reach of the e + e~ collider, the difference between the SM 
and MSSM cross sections is very small and will be hard to be detected. For the heavy 
CP-even and the CP-odd Higgs bosons, the cross sections are rather small especially for 
large values of tan (3. Only for small values of tan f3 and relatively small A and H masses 
that the cross sections exceed 10~ 2 fb in the unpolarized case. 
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APPENDIX: Scalar Functions 



Our conventions for the scalar Passarino-Veltman functions and for the tensor integrals 
are as follows \p is the renormalization scale and D the space-time dimension] 



/ , 2\ 4 j~) f d k 1 



16tt 5 



A(m 2 ) = /i 4 



(2tt) d k 2 - m 2 



(A.l) 



16tt 2 



B {plmlml)=^- D f 



d u k 



(2tt) d [k 2 - m\][{k + Pl ) 2 - ml] 



(A.2) 



16vr 2 



ri / 2 2 2 2 2 2\ 



A-D 



d D k 



1 ■ h ■ h h 



(2n) D [k 2 - ml] [(k + Pl ) 2 - ml] [(k + Pl + p 2 ) 2 - ml] 



(A.3) 



16n'- 



7-, / 2 2 2 2 / , \2 / , \2 2 2 2 2\ 

:D 0; ^ u {p 1 ,p 2 ,p 3 ,p 4 , (pi + p 2 ) AP2+P3) ;m ll m 2 ,m 3 ,m 4 ) = /J 



4-D 



d u k 



1 ■ h ■ h h 



(2-k) d [k 2 - m\] [(k + pi) 2 - ml] [(k + p 1 + p 2 ) 2 - mj] \(k + Pi + p 2 + P3) 2 - mj] 



(A.4) 



with the tensor decomposition 

C" = pZCu+pZCu 

= P^C 2l +p^ 2 C 22 + {p^ 2 +p^ 2 )C 23 + g^C 24 
D» = rfD n +p%D 12 +p%D 13 
£>"" = p1plD 21 +p^p v 2 D 22 +p^D 23 

+(pfp5 + p^)D 24 + (pfpg + p^ 3 )D 25 + (ptp" + p^ 3 )D 26 + gr D „ 



(A.5) 
(A.6) 
(A.7) 



(A.8) 



The analytical expressions of the scalar functions can be found e.g. in Ref. | 28fl . 
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Table 1: 





h 


H 


A 


g$w 


sin(/3 — a) 


cos(/3 — a) 





g®Az 


cos(/3 — a) 


— sin(/3 — a) 





g<s>H ± w ± 


=F cos(/3 — a) 


± sin(/3 — a) 


1 


g<$>G+G- 


cos(2/3) sin(/3 + a) 


cos(2/3) cos(/? + a) 






Tab. la: The Higgs-vector boson couplings g$vv [normalized to the SM Higgs coupling 



9h°vv 



V2d 



1/2 



My], and the Higgs-Higgs-vector boson couplings [normalized to gw 



(\/2Gf) 1 I 2 Mw and gz = {V^Gf) 1 ^ 2 Mz for the charged/neutral weak couplings]; the latter 
come with the sum of the Higgs momenta entering and leaving the vertices. 





g$uu 


g^dd 


h 


cos a/ sin (3 


— sin a/ cos (3 


H 


sin a / sin f3 


cos a 1 cos (3 


A 


1 / tan f3 


tan (3 



Tab. lb: Higgs boson couplings in the MSSM to fermions and gauge bosons relative to the SM 
Higgs couplings. 



g L ' R /$ 


h 


H 




„R 

y $>xtxj 


Qjj sin a — cos a 
Qij sin a — S^ cos a 


—Qji cos a — Sj { sin a 
—Qij cos a — sin a 


-Q* { sin(3 - S* iC os(3 
Qij sin /3 + cos /3 




Q'j* sin a + S'-* cos a 
Q'lj sin a + SL cos a 


—Q"j* cos a + Sj* sin a 
—Q"ij cos a + S'-j sin a 


-Qy* sin/5 + S^* cos/5 
sin (3 — S'-j cos /3 



Tab. lc: The couplings of the neutral Higgs bosons to charginos and neutralinos, normalized 

r _ il/2 „ „ _ 

to g\y = ylGp Mw\ the matrix elements Qij /Sij and Q^/S^ can be found in Ref. Iftl]] . 
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Table 2: 
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$ 
9i 


$ 

92 


#3 






u 
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cos aj sin/? 


— sm(ct + p) 


— sin a / sin (3 


cos q, / sin /3 


u 


H 


sin a / sin (3 


cos(o; + /3) 


cos a/ sin/3 


sin a / sin p 1 




A 








1 


-l/tan/3 




h 


— sin a/ cos (3 


— sin(o; + (3) 


cos a / cos (3 


— sin a/ cos p 1 


d 


H 


cos a 1 cos /5 


cos(o; + /3) 


sin a / cos (3 


cos a/ cos /3 




A 








1 


— tan (3 



Tab. 2a: Coefficients in the couplings of neutral Higgs bosons to sfermion pairs. 











Ijlj 


Shlilj 


9m l l j 


9 Aki j 




(2s 2 w - 1) sm(P + a) 


-(2s 2 w - l)cos((3 + a) 





&r&r 


2s 2 w sin(/3 + a) 


-Isly cos(P + a) 





VhVL 


sin(/3 + a) 


— cos(/3 + a) 






Tab. 2b: The couplings of the neutral Higgs bosons to left- and right-handed sleptons in the 

1/2 n 



absence of mixing, normalized to g' w = \/2Gf 



w 
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Figure Captions 



Fig. 1: Feynman diagrams contributing to the process e + e~ — > in the Standard Model 
(a), the additional diagrams for the production of the MSSM CP-even Higgs bosons 
e + e~ — > $7 with $ = h, H (b), and of the CP-odd Higgs boson e + e~ — > A7 (c). 

Fig. 2: The total cross sections for the production of the Standard Model Higgs boson, 
e + e" — * 7i? , at two center of mass energies y/s = 500 GeV and ^/s = 1.5 TeV as 
a function of the Higgs boson mass. The solid, dashed and dotted lines are for the 
unpolarized, left-handed polarized electrons, and left-handed polarized electrons 
and right-handed polarized positrons, respectively. 

Fig. 3: (a) The total production cross section <r(e + e~ — > 'yH ) for unpolarized beams as 
a function of for M H = 100 (solid), 150 (dot-dashed), 200 (dashed) and 300 
(dotted) GeV. (b) The differential cross section da/dcos9 at a center of mass energy 
y/s = 500 GeV for M H = 100 GeV for unpolarized and polarized e + e~ beams. 

Fig. 4: The total cross section for the production of the lightest MSSM Higgs boson, 
(j(e + e~ — > 7/1), in the unpolarized case at y/s = 500 GeV as a function of M h 
for tan/3 = 2 and 50. The solid lines include all MSSM contributions, while the 
dashed lines include only the contributions from the standard and Higgs sectors 
(2HDM). For the SUSY parameters we have: in (a) M 2 = 200 GeV, /i = 200 GeV, 
Mi = 200 GeV, M~ q = 500 GeV, and in (b) M 2 = 500 GeV, n = 500 GeV, M~ x = 500 
GeV, Mq = 500 GeV. 

Fig. 5: The total production cross section cr(e + e~ — > 7/1) for unpolarized beams at y/s = 500 
GeV as a function of tan (3 for M A = 50 GeV and M A = 500 GeV. The scenarios 
for the MSSM parameter space are as in Fig. 4. 

Fig. 6: The difference between the MSSM and the SM cross sections at cos 8 = normalized 
to the SM values, as a function of the SUSY parameters M 2 and /i for tan (3 = 2 
and 50. For the other parameters, we take M A = 1 TeV and Mg = Mj = 500 GeV. 

Fig. 7: The total production cross section a(e + e~ — > 7/1) as a function of the cm. energy 
v/i for M A = 50 GeV and M A 500 GeV and tan/? = 2 and 50 in the MSSM (solid 
curves) and in the 2HDM (dashed curves). The SUSY parameters are chosen as: 
M 2 = 500 GeV, fi = 500 GeV, M l = 500 GeV and M~ q = 500 GeV. 

Fig. 8: The total cross section for the production of the heavy CP-even Higgs boson in the 
MSSM, <r(e + e~ — > ^H), in the unpolarized case at y/s = 500 GeV as a function of 
M H for tan/5 = 2 and 50 in the MSSM (solid curves) and 2HDM (dashed curves). 
The SUSY parameters are as in Fig. 4. 
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Fig. 9: The total cross section for the production of the CP-odd Higgs boson in the MSSM, 
<x(e + e~ — > 7A), in the unpolarized case, at ^/s = 500 GeV as a function of Ma for 
tan/3 = 2 and 50 in the MSSM (solid) and 2HDM (dashed). The SUSY parameters 
are as in Fig. 4. 
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Fig. 7 
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